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Abstract 

The physical principles and the mathematical structure involved in deriving an 
analytical representation of the internal structure of the sun are discussed. For 
a two-parameter family of a non-linear matter density distribution, the run of 
mass, pressure, temperature, and luminosity throughout the sun are expressed in 
terms of Gauss ' hypergeometric function. The system of differential equations 
governing hydrostatic equilibrium and energy conservation for the sun generates 
another field of application of special functions. 

1 Introduction 

The structure of the sun is determined by conditions associated with mass 
conservation, momentum conservation, energy conservation, and specific modes 
of energy transport through the sun. When considering its internal structure, 



rotation and magnetic fields can be neglected since the sun is spherically sym- 
metric. It may come as a surprise that much of what has been observed and 
theorized about the sun can be accounted for in terms of very basic physical 
laws: Newton's laws of gravity and motion, the first two laws of thermodynam- 
ics, Einstein's law of the equivalence of mass and energy, Boyle's law, Charles' 
law of perfect gases, and Heisenberg's uncertainty principle. The outline of the 
theory of the structure and evolution of the sun has been formulated in the 
first half of the twentieth century; it is connected with names like Lane, Em- 
den, Schwarzschild, Eddington, Chandrasekhar, Hoyle, and Fowler, see Mathai 
and Haubold (1988). In the second half of this century, theories about the sun 
were greatly refined; this is due in part to new observation techniques and to 
computer simulations of its structure and evolution. By and large, the overall 
picture of the structure and evolution of the sun seem to be well understood. 
However, there is a serious discrepancy between the theory of how the sun shines 
and the most direct experimental test of this theory. This discrepancy is more 
specifically called 'the solar neutrino problem'; it refers to the fact that the sun 
is a volume source of neutrinos which are particles produced by thermonuclear 
reactions in the deep interior of the sun, and that the copious flux of solar neu- 
trinos predicted by theory does not match the flux detected by experiments on 
earth over the past 25 years, see Abdurashitov et al. (1994), Anselmann et al. 
(1994), Davis (1993), and Nakamura (1993). This problem has been studied by 
Davis and Bahcall among others, see Bahcall (1989). Solar neutrinos are the 



only particles that have the ability to travel from the center of the sun to its 
surface almost without interaction with solar matter, escaping freely into space 
carrying the most direct information about physical processes in the deep solar 
interior. The solar neutrino problem, which shall not be the subject of this 
paper, stimulated further studies on solar models, both by employing modern 
computing tools and using the analytical techniques of mathematics. This was 
also the justification for an effort to reconsider the derivation of analytic so- 
lutions to the system of differential equations of solar structure based on the 
very basic physical laws mentioned above. Particularly, the solar neutrino prob- 
lem has been considered to be a reason to pursue more actively the problem of 
obtaining analytic formulae which give a description of the gravitationally sta- 
bilized solar fusion reactor, thus showing that methods of the integration theory 
of generalized special functions applied to solar physics constitute another field 
of application for these functions, see Mathai and Haubold (1988). 

The nuclear reactions which cause the sun to evolve are sufficiently slow that 
the sun may be assumed to pass through a series of equilibrium configurations. 
The model for the current internal structure of the sun may be thought of as 
representing the sun at an instant of time. Separating the time dependence 
of the evolution of the sun from the equations governing its internal structure 
allows one to replace the time dependent partial differential equations by four 
simultaneous, non-linear, ordinary differential equations of the first order. The 
four equations represent the radial gradients of mass, M(r), pressure, P(r), 



temperature T(r), and luminosity, L(r), where r denotes the radial distance 
from the centre of the sun. Since there are four equations but more than four 
unknown physical variables, one needs additional constitutive equations, before 
the system can be solved: an equation of state for solar matter, a nuclear energy 
production rate, and an opacity law. The full system of equations must be solved 
subject to at most four boundary conditions at the surface and the centre of the 
sun. These boundary conditions ensure that the structure of the sun can be cal- 
culated from the four differential equations but they do not ensure that there is 
a single unique solution, see Chandrasekhar (1939) and Stein (1966). The proce- 
dure of numerically integrating the solar structure equations takes advantage of 
large electronic computers and makes it possible to include a variety of detailed 
physical effects and to vary parameters at will, see Bahcall (1989) and Noels et 
al. (1993). A second procedure providing a solar model, whose contents can be 
understood intuitively without resorting to numerical techniques, is based on 
Buckingham's theorem which states that a system characterized by n physical 
variables can be described by a set ensemble of n-r dimcnsionslcss products of 
variables, where r is the number of variables whose dimensionless representations 
are linearly independent. This approach provides a qualitative explanation of 
the fundamental stellar structure equations through dimensional analysis. This 
analysis suggests that more detailed physical and mathematical theories are es- 
sentially theories of factors of proportionality. They eventually yield numerical 
values for these proportional factors because more physical assumptions have 



to be made, see Bhaskar and Nigam (1991). The third procedure treats the 
solar structure equations by rigorous mathematics leading to the Lane-Emden 
equation which is a second-order non-linear differential equation describing the 
structure of a polytrope gas sphere. However, explicit analytic solutions of the 
Lane-Emden equation exist only for the values n=0,l, and 5 of the polytrope 
of index n, not covering the specific physical model for the internal structure of 
the sun, see Chandrasekhar (1939) and Horedt (1990). 

In addition to those three procedures for constructing solar models, there 
exists another approach for finding solutions to the solar structure equations 
which consists of making a specific assumption for straightforward analytical 
integration of these equations. It is possible to obtain analytic solar models 
by separating the hydrostatic component from the energy-transport component 
of the structure equations. For that purpose, an analytic density distribution 
assumption, namely, that the matter density in the sun varies non-linearly from 
the center to the thought surface, where the density goes to zero, must be made 
see Stein (1966) and Mathai and Haubold (1988). Itis then possible to inte- 
grate the equations of mass conservation, hydrostatic equilibrium, and energy 
conservation through the sun. Together with the equation of state of a perfect 
gas, the run of density, p(r), mass, M(r), pressure, P(r), temperature, T(r), 
and luminosity, L(r), are determined and can be derived in the form of analytic 
formulae. The physics of the problem requires only three independent boundary 
conditions: M(r) — > and L(r) — > at radial distance r = 0; T — > Tq = and 



p — > po = at the radius r = Rq of the gaseous configuration. The requirement 
that p and T tend simultaneously to specific values, in this case zero, is only one 
condition since the point at which this occurs is arbitrary. This ambiquity can 
be removed, in principle, by assigning the total mass. The boundary is required 
to be at the point where M(r = Rq) = M© and this provides the fourth condi- 
tion. Hence, the central density, pressure, temperature, and total rate of energy 
generation are determined as a function of the sun's mass and radius. How- 
ever, by assuming an analytic matter density distribution, the energy-transport 
equation of the system of structure differential equations can be satisfied at only 
one typical point in the sun. The procedure thus established to construct an 
analytic model of the solar interior allows to determine the factors of propor- 
tionality which remain to be an open problem in the dimensional analysis. The 
procedure also reveals that the run of all physical variables for the solar model 
can be expressed in terms of Gauss' hypergeometric function. 



2 Matter Density Distribution 

For the integration of the system of differential equations governing the in- 
ternal structure of the sun one has to make a choice for an unknown function 
that still leaves room for physical justification of this choice. By intuition one 
expects that mass is an increasing function while density, pressure and temper- 
ature are decreasing functions from the centre to the surface of the sun. Thus 



we make the working hypothesis that the matter density distribution p(r) varies 
with the distance variable r, as 

Si 7 



(2.1) 



p(r) = Pc 
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Rq 



,6>0, 7 >0, 0< — <1, 

Hq 



where S and 7 are kept as free parameters to ensure that the density distribution 
determines properly the mass, pressure, and temperature distributions in the 
sun. Equation (2.1) takes into account that the chosen density distribution 
reflects the central value of the density p(r = 0) = p c and satisfies the boundary 
conditon p(r = Rq) = 0, where Rq denotes the solar radius. Also, equation 
(2.1) implies that p oc M & /R & , where the constant of proportionality depends 
only on the radial mass distribution and the radial distance. 



3 Distribution of Mass 

If M (r) represents the total mass contained within the radius r, and p(r) is 
the density at r, then 

(3.1) ^W =4H r»p(r). 

dr 

Using the assumed non- linear density distribution in equation (2.1), integration 
of (3.1) throughout the Sun leads to the distribution of mass 

where 2^1 (•) denotes Gauss' hypergeometric function, which contains the pa- 
rameters 5 and 7 of the matter density distribution in equation (2.1), see Luke 



(1969) and Mathai (1993). Equation (3.2) satisfies the boundary condition 
M(r = 0) = and can be used to determine the central value of the matter 
density distribution p c in equation (2.1) in terms of the parameters 8 and 7 of 
the chosen model of the sun. The condition M(r = i? Q ) = M Q in equation 
(3.2) reveals that 

(3 ' 3) Pe = ^7iUA* )' 

if 7 in equation (2.1) is kept as a positive integer; (3.3) is obtained by using 
the following relation for Gauss' hypergeometric function of argument one, see 
Luke (1969) and Mathai (1993): 

2 F 1 (a, b; c; 1) = T(c)T(c -a- b)/T(c - a)T(c - b). 

Equation (3.3) can be used to select the appropriate values of the parameters 
8 and 7 specifying the solar model with matter density distribution given in 
equation (2.1). 



4 Distribution of Pressure 

If g = GM(r)/r 2 is the gravitational force per unit mass at r due to the 
attraction of the mass interior to r, then 

dP(r) GM(r)p(r) 
{ ' dr r 2 

is the equation of hydrostatic equilibrium of the spherical self- gravitating sun 
with dP{r)/dr being the pressure gradient. The internal pressure produced by 



the weight of the overlying layers increases towards the centre while the gas 
and radiation pressure must increase correspondingly to achieve the balance of 
forces for equilibrium. This increase is obtained by inward increases of both 
temperature and density. The internal pressure needed to achieve a balance is 
the gravitational force per unit mass (GMq/Rq) times the mass per unit area 
(Mq/Rq) which gives that P oc G Mq/Rq for any spherical body in hydrostatic 
self- gravitation, where the constant of proportionality is again determined by 
the radial distribution of mass in the sun and the particular radial distance at 
which P is measured. The constant of proportionality can be determined by 
integrating equation (4.1) throughout the volume of the sun by using equations 
(2.1) and (3.2) for the density and mass distribution, respectively. We obtain 



where 2-F\(-) is Gauss' hypergeometric function and (— -f) m = T(— 7+m)/r(— 7) 
is Pochhammer's symbol that often appears in scries expansions for hypergeo- 
metric functions. The solution of equation (4.1) given in equation (4.2) confirms 
the condition P(r = Rq) = and gives the central value of the pressure accord- 
ing to the chosen solar model characterized by S and 7 in equation (2.1): 




(4.2) 
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5 Temperature Distribution 

The simplest theory of solar structure is that of a polytrope. These solar 
models obey an equation of state of the form P = ifp( n+1 )/™ throughout the gas 
sphere. Since temperature does not explicitly appear in this relation between p 
and P, equations (3.1) and (4.1) may be solved independently of the tempera- 
ture and luminosity gradients. This equation of state leads to the Lane-Emden 
equation for polytropic gas spheres, which is an ordinary differential equation of 
second order, but can be reduced, by suitable transformations of the variables, 
to an equation of the first order (Chandrasekhar, 1939; Horedt, 1990). In the 
sun the density is so low that at the temperatures involved the solar material 
behaves almost as a perfect gas, having molecular weight p and obeying the 
perfect gas law 



where k is Boltzmann's constant and Na Avogadro's number. Substituting 
p oc Mq/Rq and P oc GMq/Rq in (5.1) reveals the dependence of temperature 
on the mass and radius of the sun T oc yuM Q /_R Q , where the constant of propor- 
tionality depends on the mass distribution and the radial distance. We obtain 



the detailed temperature distribution throughout the sun by using equations 
(2.1) and (4.2) to rewrite the equation of state given in (5.1), that is: 



Equation (5.2) satisfies the boundary condition T(r — R Q ) = and allows to 
determine the central value of temperature of the sun as a function of 5 and 7 
contained in equation (2.1): 



The procedure is to construct a model for the internal structure of the sun by 
assuming a matter density distribution and subsequently integrating the system 
of differential equations. At this point two remarks are in order. The contribu- 
tions of the radiation pressure to the total pressure and the radial dependence 
of the mean molecular wcigth /i have been neglected in equation (5.1). 

The total pressure P at any point in the sun is the sum of the gas pressure 
and the radiation pressure, P = P g + P r , where P g is given in equation (5.1) 
and P r = ^aT 4 , where a is a constant. Writing P g — f3 P and hence P r = 
(1 - (3)P, it follows that P = aT 4 /3(l — (3). This ratio of radiation pressure 
to gas pressure increases towards the center of the Sun but even there the gas 




(5.2) 





(5.3) 




pressure exceeds the radiation pressure by more than three orders of magnitude. 
This justifies neglectingt the radiation pressure has been neglected in equation 
(5.1) (Chandrasekhar, 1939). 

The outward flow of energy inside the sun is driven by the temperature 
gradient and resisted by the opacity of the material. The nuclear energy gen- 
erated within the sun has to continually replenish that radiated away from the 
surface. This energy generation by nuclear reactions causes the solar chemical 
composition to change and keeps the sun evolving. Since the gas inside the sun 
is completely ionised, the mean molecular weight /j, in equation (5.1) is given 
by \i = (2X + jY + \Z)~ X , where X, Y, Z are relative abundances by mass of 
hydrogen, helium, and heavy elements (X + Y + Z = 1). The dependence of 
X, Y, Z on the radial distance variable r, which is governed by kinetic equations, 
cannot be determined in the procedure of constructing an analytic solar model 
by assuming a matter density distribution as given in equation (2.1). Thus, the 
mean molecular weigth has to be treated as constant in the following. This as- 
sumption does not reflect the situation in the real sun because nuclear reactions 
have changed the originally uniform chemical composition throughout the sun. 



6 Nuclear Energy Generation Rate 

Nuclear energy production in the sun depends heavily on the temperature 
of the material and is very concentrated towards the center of the sun. This ex- 



plains that calculations of the internal structure of solar-type-stars made consid- 
erable progress even before the physical mechanism of the production of energy 
by nuclear reactions was understood. The rate of nuclear energy generation can 
be written (Mathai and Haubold (1988), as 



where eo is a physical constant depending only on the chemical composition of 
the solar material and the units chosen. Substituting p(r) and T(r) in equation 
(6.1) by equations (2.1) and (5.2), respectively, and taking advantage of the 
fact that < (^— ) < 1, the energy generation rate in equation (6.1) can be 
represented in the form of a polynomial 



(6.1) 



e(p,T) = e p n (r)T m (r) 



(6.2) 




where / denotes the expression 
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no + n\ + . . . + n 2l = q 



and 



(6.4) 




Note that the representation of e(r/Ro) in equation (6.2) is essentially deter- 
mined by the four free parameters 5, 7, n and m in equation (2.1) and (6.1), the 
only restriction being that 7 be a positive integer. The coefficients ao, a\, . . . a2 7 
in equation (6.3) are determined by the following polynomial of degree 27 in 

, <5 
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7 Luminosity Function 



Let L(r) be the function representing the flow of integrated radiation across 
a sphere of radius r. If e is the energy produced per unit time by nuclear 
reactions in each unit mass of solar material, then the balance between energy 
generation inside the sun and energy loss through its surface is governed by the 
equation of energy conservation 

dL(r) 



(7.1) 



dr 



4nr p(r)e(r), 



where e(r) is given by equation (6.2). Integrating equation (7.1) over the sun's 
interior leads to the luminosity function in terms of Gauss' hypergeometric func- 
tion: 



= A^p n c +1 T™R% 
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s* = s + -(3 + 2gr) + m + 2n 2 + . . . + (27)n 27 , 



where / is given in equation (6.3) with s substituted by s* . Equation (7.2) 
satisfies the condition L(r = 0) =0 and gives for the total energy output 
L(r = Rq) = L Q , 

(7.3) L Q =Ave oP ^T™R%\f- 
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8 Conclusions 

Assuming a two-parameter family of matter density distributions in equa- 
tion (2.1) made possible the analytic integration of the differential equations for 
conservation of mass, momentum, and energy throughout the sun, respectively 
equations (3.1), (4.1) and (7.1). This procedure shows that hydrostatic equi- 
librium, equation of state, and overall energy conservation determine the state 
of the central solar conditions, particularly the gravitationally stabilized solar 
fusion reactor. The mathematical method chosen reveals the factors of pro- 
portionality which are kept undetermined in dimensional analysis, commonly 
pursued to understand astrophysical relationships between global parameters of 
the sun. A common mathematical element of the derived distributions of mass 
(3.2), pressure (4.2), temperature (5.2), and luminosity (7.2) throughout the sun 
is Gauss' hypergeometric function 2-Fi(-) which is numerically easily accessible 
through mathematical programs like Mathematica (Wolfram, 1993). 



It has been emphasized above that the assumption of an analytic matter 
density distribution means that the equation for the transport of energy by 
radiation through the sun can be satisfied at only one specific point in the sun. 
The flow of radiant energy per unit area through the sun is proportional to the 
ratio of radiation pressure gradient and opacity per unit volume. That is 

d^aT^/dr THTjdv 

(8.1) Hoc — 2 oc , 

Kp Kp 

where k denotes the opacity per mass unit at temperature T and density p. 
Because the energy flowing out through the sun is transported by radiation, we 
find for the luminosity L & 

, s 2 RlT 3 dT/dr 

(8.2) L Q oc AttR%H oc -2 '- — . 

Kp 

Since for a given solar structure p oc M Q /R Q and T oc Mq/Rq, it follows for 
Lq that 

(8.3) Lq oc -M|. 

For solar the composition, Kramer's power law approximation for the opacity 
given by 

(8.4) k oc k qP T- 7/2 , 

where «o is a physical constant depending on the chemical composition of the 
solar material and the units chosen, which leads to a luminosity-mass-radius 
relation, 

(8.5) Lq oc Af" /2 i?Q 1/2 . 



The differential equation governing the outward flow of energy driven by the 
temperature gradient and resisted by opacity in (8.1) is 



where c denotes the velocity of light. To satisfy the radiative energy transport, 
taking into account the density distribution assumed in equation (2.1), the tem- 
perature gradient dT/dr in equation (8.6) has to be equal to the temperature 
gradient in equation (5.2). This condition can be satisfied at only one specific 
point in the solar interior, for example at the boundary of the nuclear energy 
producing core region (at r « 0.3i? Q where L « L & ). 
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